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Abstract 

In the present paper our aim is to introduce some models for the generalization of the kinetic 
theory of electrons and phonons (KTEP), as well as to study equilibrium solutions and their 
stability for the generalized KTEP (GKTEP) equations. We consider a couple of models, relevant 
to non standard quantum statistics, which give rise to inverse power law decays of the distribution 
function with respect to energy. In the case of electrons in a phonon background, equilibrium and 
stability are investigated by means of Lyapounov theory. Connections with thermodynamics are 
pointed out. 
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INTRODUCTION 



As pointed out by Koponen, fractal or inverse power law distributions are of interest in 
modeling various meaningful situations in solid state physics [lj. An example, treated in 
ref . |l( , is the thermalization of a nonequilibrium electron-phonon system. 
Until recently, however, there has been little guidance on how to generalize the kinetic the- 
ory of electrons and phonons obeying non-Gibbsian statistics. In a very recent paper Q] a 
generalized kinetic theory for electrons and phonons has been proposed. Modified collision 
terms were introduced, in order to allow applications not only to electrons, but also to other 
particles (obeying a general statistics) which interact with a crystal lattice. 
In the present paper our aim is to introduce some models for the generalization of the KTEP 
equations as well as to study equilibrium and stability for the GKTEP equations. We in- 
troduce a couple of models for the GKTEP, relevant to non standard quantum statistics 
recently proposed in the literature, and compare them from the point of view of Kinetic 
Theory. Both of them give rise to inverse power decays of the distribution function with 
respect to energy. In absence of external disturbances, it is possible to treat phonons as 
background at equilibrium. This assumption implies a drastic simplification of the GK- 
TEP, since we have just to deal with the equation for electrons. Equilibrium solutions and 
their stability are investigated by means of Lyapunov's theory. The connection between 
mathematical results and Thermodynamics is shown. 

First we recall the GKTEP equations and their stationary solutions. Let N g = N g (k, x, t) 
be the distribution function of phonons (quasi- momentum k, energy uo g {k) of branch g of 
the phonon spectrum and n p = n(p, x, t) the distribution function of electrons [quasi- 
momentum p, energy e(p)]. 

The generalized kinetic equations for phonons and electrons read Q]: 
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where the w's are transition probabilities and electron-electron interactions have been ne- 
glected. Observe that in the limit 

*( Ng ) _> i + i v 9 , $(iV s ) - iV s , (6) 

^(n p ) -> 1 - n p , <£>(n p ) -> n p , (7) 

the Bloch-Boltzmann-Peierls (BBP) equations Q are recovered. The non negative functions 
$ for phonons [N g G [0, oo)], and if), (p for electrons (n p G [0, 1]) represent the arrival 
and departure state availability, respectively. Moreover we assume: 
i) $>/ty and tp/ip are monotonically increasing functions of N g and n p , respectively Q]; 
m) the following conditions hold 

V(0)=0, $(0)=0, (8) 
^(0) = 1, *(0) = 1, (9) 
V>(1) = 0, ^(+oo) = +oo. (10) 

Conditions (jHJ) mean that no transition occurs when the initial state is empty, conditions 
(JHJ) mean that the emptiness of the arrival state neither inhibits nor enhances a transition 
and conditions (jlUj) are characteristic of fermions and bosons, respectively. 
In 3] it has been shown that, implicitly, the equilibrium solutions of the generalized system 
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ln#| = -^, (12) 

(* means "at equilibrium"), where T is the absolute temperature of the mixture and \x is 
the chemical potential of the electron gas. 



MODELS FOR THE GENERALIZATION 

Preliminarily, consider certain particles, whose equilibrium distribution function is /, 
which obey a statistics described by the parameter A G [—1, +1]. Hereinafter hatted quan- 
tities are referred to these particles. In particular, for A = —1, 0, +1, we have fermions, 
classical particles, and bosons, respectively. Suppose, we invert this distribution function 
with respect to exp(— £) [where £ = (e — fi)/T] than we obtain: 

exp(-£ ) = F(f, A) . (13) 

In particular, for electrons £ = (e p — /i)/T, for phonons £ = uj 9 /T, since, in this case, ft = 0. 
The statistics of the particles we deal with is described also by the couple of function 0(f, A) 
and ip(f, A), essentially positive, and obeying the following conditions: 

0(0, A) = , ^(0, A) = 1 . (14) 

In particular we have 

<p(Pv) = 0( n P> - 1 ) , 4>(Pv) = ^( n pi _1 ) > ( 15 ) 

9(N g ) = 0(N g , +1) , 9(N g ) = ${N g , +1) . (16) 

From (jllj) and (|12jl we know that at equilibrium 

j = exp(-0 . (17) 



Now, Eqs. (USD and JTZJ) give 



7 = A) . (18) 



Let us introduce the following assumptions: 

i) is given by Eq. (j!8|) not only at equilibrium; 

ii) 0(f, A) does not depend on A, and we take <p(f, A) = 0(f, 0); 
Hi) ^(/, 0) = 1. 

Properties i), ii) and Hi) actually hold for standard statistics (Fermi-Dirac, Maxwell- 
Boltzmann, and Bose-Einstein) and, more in general, we postulate them. From Eq. (|18|). 
by accounting for ii) and Hi), we get 

= F(f, 0) , 4, = ^| . (19) 

We introduce now two models which account for an inverse power law decay of the 
distribution function with respect to energy: 



1) The first model that we consider has been proposed by Buyiikkilig et al. |3, I^J . Relevant 



applications to the blackbody problem can be found 
g-deformed exponential and logarithm []]]: 
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Preliminarily, we define the 



exp q (x) = — — , \n q (x) = L/J— . 20 

[l-(q-l)x]— 9-1 

The deformed quantum distribution function is given by 

/ = tt — r » (21) 

l/exp ff (-£) - A 

where the parameter q is greater or equal to 1 (in the last case the standard statistics are 
recovered). Since F(f, A) = exp{ln 9 [l/(l// + A)]}, from Eq. (|T9*j) we get 

<p = exp(ln (? /) , tp = exp |hi g / - \n q - ^ X ^ J ■ ( 22 ) 

We have limj^ y? = and limj^oip = 1 for q < 2, otherwise it either vanishes (for A < 0) 
or it diverges (for A > 0), 

^(1, A) = exp(-ln g — j— ) , (23) 



and 



1 + A 



lim ip(f, A) = exp ( ) , for A > . (24) 
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FIG. 1: a) and b): plot of <p and tfj vs n and q for the model of Biiyiikkilig et al; c) and d): plot of 
ip and ip vs n and k for the model of Kaniadakis. 



For fermions we observe that ^(1) = is not fulfilled. Moreover, in the case of bosons the 
condition ^(+00) = +00 does not hold. 

Figures la and lb describe the couple of function ip and ip versus both n and q. Figures 
2a and 2b describe the couple of functions $ and \1/ versus both N and q. Figure lb shows 
clearly the violation of condition ip(l) = 0. Both ip and $ start with a vanishing derivative, 
while both ib and \1/ start with an infinite derivative. 

2) Recently K aIli a dakis Q Q has proposed a new ela.ica, distribution function, together 
with its quantum counterparts. An application to Bose-Einstein condensation can be found 
in ^(|. Preliminarily, we define the new K-deformed exponential and logarithm: 
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Observe the remarkable property exp| K |(— x) = 1/ exp| K }(x). The deformed quantum dis- 
tribution function is given by 

1 



(26) 



exp w (£) - A 

where the parameter k is taken here greater or equal to (in the last case the standard 
statistics are recovered). Since F(f, A) = exp[— hi{ K }(l// + A)], from Eq. (fTTJJ) we obtain 

1 + AA 



tp = exp(ln w /) 



ijj = exp ln {K} / + ln {/t} 



./ 



(27) 



We have \imf_^ o = and lirn/^oV' — 1 f° r K < 1> otherwise, it either vanishes (for A < 0) 
or it diverges (for A > 0), and 



^(1, A)=exp[ln w (l + A)] . 
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FIG. 2: a) and b): plot of $ and f vs JV and q for the model of Biiyiikkilig et al; c) and d): plot 
of <E> and vs N and k for the model of Kaniadakis. 
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For electrons, the last result agrees with the condition if>(l) = 0. In the case of bosons, the 
condition for $ when x — > oo is fulfilled as well. Figure lc and Id describe the couple of 
functions tp and if) versus both n and k. Figures 2c and 2d describe the couple of functions $ 
and \& versus both N and k. Figure Id shows clearly the fulfillment of conditions ^(l) = 0. 
Both if and $ start with a vanishing derivatives, while both if) and \I/ start with an infinite 
derivative. 



ELECTRONS IN A PHONON BACKGROUND: EQUILIBRIUM AND STABILITY 

A usual approximation to the BBP equations consists in considering the relaxation due to 
p-p interactions to be much quicker than the one due to e-p interactions. This assumption 
amounts to fix, in the kinetic equation for electrons, the distribution of phonons as the 
equilibrium function (jllj) at a given temperature T. 
The kinetic equation for electrons reads now 

^E + „.^E_ eJ5 .^E = (^y , (29) 
at ox op \ ot J ep 

where o means that N* has been substituted for N g . In order to study equilibrium and its 
stability, it is useful to introduce the following functional: 



\ / ep 

which can be rewritten as follows: 



exp 



dp , (30) 



C = 2 / / w ep (p', k^p)5(e p - e p / - u g ) In 



xWN;)ip(n p ,)if;(n p )-<p(n p )iP(n pl )*(N;)}^ < 0, (31) 
since (A- B)\n(B/A) < 0. 

In the space homogeneous and forceless case, equilibrium solutions are given by 

(On 
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ep 



Proposition 1. 

Eq. Rflfy is equivalent to 



$(iV!) <pM iP(n p ) = if(n p ) iP(n p/ ) V(N* g ) , V p , k . (33) 



Proof. In fact, we observe that Eq. (|3*3*Jl =^> Eq. [see Eq. (j3J)]. On the other hand, 
from Eq. (JSOD it follows that Eq. (J32D =>• C = 0. Tanks to Eq. flUU), we obtain C = =>• 

Eq. d22D • 

Eq. (j22I)) by taking into account Eq. (fTTj) and e p = e p i + u g , shows that, at equilibrium, 
ha[(p(n p )/ip(n p )} + e p /T is a collisional invariant for electrons, that means n* is given by 
Eq. (JT2J). In order to study the stability of such an equilibrium solution, we introduce the 
following functional 

2 r 

L = H e + - J e p n p dp , (34) 

where H e = J H e (n p ) dp and &H e /dn p = 2 \rt((p/ip) (observe that H e is a convex function 
of n p , since tp/ifi has been assumed to be monotonically increasing). 
We can prove an H theorem for the present problem: 
Proposition 2. 

L is a Lyapounov functional for the present problem. 
Proof. By utilizing the definition of L together with Eq. (|29J). one obtains L = C < 0, 
where the dot means the time derivative. Moreover, since (&H e /dn p )* = 2 (/i — e p )/T by 
taking into account electron conservation, we have 

'mu 



l\~o~\ ( n P ~ n *v) d P = ~t / e pK ~ n l) d P > ( 35 ) 



so that we can write 



l-l* = in 



K+(^]\n p -n;) 



dn p 



dp > , (36) 



due to the convexity of 7Y e • 

Let us now interpret this result on a physical ground. First, we define the concentration 
N, the energy density E e , and the entropy density S e of the electron gas as follows: 

N=^/2n p dp, E e = ^-j2n p e p dp, S e = , (37) 

where the factor 2 inside these integrals accounts for degeneracy. Observe that this definition 
of ^ is consistent with the definitions of T and fi, which we have already given. In fact, at 
equilibrium, the following thermodynamical relationships |3| are recovered: 

fdS e \ /I fdS e \ 1 



dNj Ee T \dE e J N T 



(38) 



Now, from the definition of L, it is easy to realize that the meaning of L < is nothing but 
the Clausius inequality S e > E e /T. 
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